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Using a three-current model for heat, spin-up and spin-down electrons, the thermodynamics of 
irreversible processes predicts that a temperature gradient gives rise to a spin current under the 
conditions used to measure what is called the spin Seebeck effect. A diffusive current proportional to 
a gradient of the chemical potential is known in thermochemistry as a Soret effect or thermophoresis. 



I. INTRODUCTION 

Spintronics is concerned with using the spin of the elec- 
tron in novel electronic devices. The spin property adds 
a degree of freedom to the electron transport. Recently, 
it has become evident that a number of studies of spin- 
dependent transport included also heat transport and the 
term "spin caloritronic" was coinecffl While the charge, 
spin and heat currents were considered by spintronics 
founders^, much attraction to the inclusion of heat in 
spin-dependent transport st em s from the discovery of the 
so-called spin Seebeck effecW^. 

The Seebeck effect is known in the theory of non- 
equilibrium thermodynamics as a cross-effect involving 
two generalized forces, a temperature gradient and an 
electrostatic potential gradienP. Here we review another 
cross-effect, the Soret effect. It is closely related to the 
Seebeck effect and we show how it may apply to spin 
caloritronics. 

Charles Soret observed that a tube containing a mix- 
ture of two salts, where one end of the tube was main- 
tained cold and the other hot, presented different salt 
concentrations at both end^. Since then, the diffusive 
current of one substance with respect to another, driven 
by a temperature gradient, has been known as the Soret 
effect or thermophoresis. In liquids, it is usually demon- 
strated using two plates, where the top one is maintained 
hot and the bottom one cold in order to avoid convec- 
tion. After a short while, the balance between the heat 
current and the diffusion flux leads the system to a sta- 
tionary stated Thermodiffusion is more commonly ob- 
served with aerosols in which the effect is quite intuitive 
: the air molecules have on average a higher velocity on 
the higher temperature side of any particle in suspension, 
thus pushing it to the lower temperature side. 



II. THREE CURRENT MODEL 

In order to examine the Soret effect in the context of 
spin transport, we consider that the charge carriers are of 
two types, up spins and down spins, labeled with a sub- 
script (+) and (— ), respectively. Then, we use a three 
current model involving a heat current density, and the 
two charge current densities. We use the notation of a 
previous worlM A similar three-current model was re- 
cently used bjEJ The three diffusive currents are the en- 



tropy current density j s and the electric current densities 
j + and j_ of the spin up and spin down carriers. Ther- 
modynamics of irreversible processes implies that there 
are linear relations between these current densities and 
their respective generalized forces, which are the Onsager 
reciprocity relations^. These relations relate the heat 
current j s and the charge current densities j+ and j_ re- 
spectively to the gradient of temperature V T, and the 
gradients of the respective electrochemical potential of 
the charge carriers, V/i+ and V/z_ according to, 
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For convenience, we adopt the notation of Valet and 
FertP^ and express the electrochemical potentials fi± as, 

fi± = ± + qV , (2) 

where q is the charge per carrier, V the applied electro- 
static potential and the mean chemical potential. In 
our discussion, we do not consider purely chemical effects, 
so V/io = 0. We identify now the physical meaning of 
the Onsager coefficients based on phenomenological con- 
siderations. 

The Onsager coefficients C-\ and £_+ describe spin 

mixing 12 . These coefficients are important when describ- 
ing the magnetic field dependence of the thermoelectric 
powei^. Here we neglect them as they simply modify 
slightly the coefficients of our final result, i.e. 

C+- = and £_+ = . (3) 

The spin-dependent electric conductivities a± of each 
spin current are determined in the absence of a temper- 
ature gradient, i.e. VT = and chemical potential gra- 
dient, i.e. V/i± =0, which then implies that ([I]) reduces 
to, 

j± = -C ±± qVV =(7±E , (4) 

where E = -V7 is the applied electric field. Thus, 
the Onsager coefficients C±± are recast in terms of the 
spin-dependent electric conductivities respectively as, 

C++ = ?± and £__ = —. (5) 
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The spin-dependent Seebeck coefficients, e±, defined 



as. 



E± = - 



q vr 



(6) 



are compared to virtuai experiments in which one im- 
poses one spin current density to vanish, i.e. j± = 0, 
which then implies that the corresponding relations in ([!]) 
reduce to, 
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Thus, the Onsager coefficients C± s are recast in terms 
of the spin-dependent electric conductivities a± and the 
spin-dependent Seebeck coefficients e± respectively as, 



and 



(8) 



The Onsager reciprocity relations state that the co- 
efficients Cij and Cji, where i,j — {s, +,— }, describe 
reverse processes under a time reversal. Since the spin 
change sign under a time inversion, up spins correspond 
to down spins in reverse processes and vice versa. Thus 
the Onsager coefficients C s ± are related to the reverse 
Onsager coefficients C TS by, 



and 
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The thermal conductivity k is determined in the ab- 
sence of spin-depend current, i.e. j + = j = 0. Using 
the relations Q and Q, the thermal equation in 
yields, 
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Thus, the Onsager coefficient C ss is recast in terms of the 
thermal conductivity coefficient k, the spin-dependent 
electric conductivities a± and the spin-dependent See- 
beck coefficients e± respectively as, 



C ss = - + (a + + (T_)£ + £_ 



(11) 



where k' is the ratio of the thermal conductivity divided 
by the temperature. Using the analytic expressions for 
the Onsager coefficents @, ^ and Q, the 

Onsager relations ([I]) become, 
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At this point, it is physically useful to define the elec- 
tric current density j and the spin-dependent electric po- 
larisation current density j p as the sum and the difference 
of the spin-dependent current densities j+ and j_ respec- 
tively, i.e. 



J- 



and 
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(13) 



Similarly, we define the effective electric conductivity a 
and the spin-dependent polarisation conductivity a p as, 



and 



(14) 



Spin relaxation occurs near the electrodes. At dis- 
tances much greater than the diffusion length the spin 
channels reach equilibrium, which implies that V (A/i) = 
and from ([2]), it follows that V/i = V/x + = V/U_. 

In this bulk limit, we define the effective spin Seebeck 
coefficient e and the polarisation spin Seebeck coefficient 
e p respectively as, 
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£ P = - q VT f ° r Jp = °- 
(15) 



The Onsager relations ( 12 ) imply for j = and j p = 
respectively that, 



(cr + e + +a-£-) VT + 
(<J+e+ - as-) V T + 
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By comparing the relations ( 15 ) and( 16 ), the analytical 
expressions for the spin Seebeck coefficients e and e p are 
respectively found to be, 
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Finally, the Onsager matrix (12) can be recast as 
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In chemistry, the Soret effect refers to the contribu- 
tion to the diffusive current of the solute with respect of 
solvent^, which is due to the temperature gradient. By 
analogy, the spin Soret effect refers to the contribution to 
polarisation current due to the temperature gradient 
VT. 

By defining the spin Soret coefficient S as, 
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it follows from (18), that the spin Soret coefficient S is 
given by, 

£ = (TpEp . (20) 



III. SPIN SEEBECK EFFECT 

We analyse now the spin Seebeck experiment. The 
platinum electrodes are assumed to produce a very strong 
spin relaxation, so that the experiment can be modeled 
with A/i = at both electrodes. Then, the diffusion 
equation for A/x implies V (A/i) = 0. As in any Seebeck 
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measurement, there is no electric current, i.e. j = 0. 



Under these conditions, using (15), the spin polarisation 



current density j p in ( 18 ) is found to be proportional to 
the temperature gradient and reduces to the elegant and 
simple expression, 
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The spin-dependent electric conductivities a± and the 
spin-dependent spin Seebeck coefficients e± can be writ- 
ten as, 



a ± = -(l±P) 



and e± = e (1 ± 77) . (22) 
Thus, to first order in /3 and 77, 



and 



s p = e I 1 



(23) 



and the expression (21) for the spin polarisation current 
density j p reduces to, 
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IrP, a spin polarisation current density j p is detected 
by the inverse spin Hall effect occurring in the platinum 
electrodes. Uchida et alFSl deduced a spin current from 
a two-current model by assuming that the temperature 
derivative of the chemical potential was spin-dependent, 
which is quite different from the present description. 

There are many systems where two types of carriers 
(of charge as well as of entropy) are invoked. For ex- 
ample, in reference to the observation of a spin Seebeck 



effect in insulating materials^, it has been conjectured 
that the difference between the effective temperature of 
magnons and that of phonons drives spin pumping into 
the electrons^. From a mere thermodynamic standpoint, 
a drift of the magnonf^ toward the cold side of the sam- 
ple is a form of thermophoresis. Phonon drag can also 
be thought of as a cross-effect between phonons and elec- 
trons^. It has been suggested that it may also play a role 
in the spin Seebeck effedP. 

As always with the thermodynamics of irreversible 
processes, the Onsager matrix expresses the existence 
of relationships among the gradients of the extensive 
state variables and their associated currents, however this 
phenomenological approach cannot address the micro- 
mechanisms responsible for the values of the matrix ele- 
ments. 



IV. CONCLUSION 



In summary, we consider the spin Seebeck effect as an 
experiment characterised by zero effective electric current 
and the equilibrium of the chemical potential at the elec- 
trodes. We show that a three-current model of heat, spin- 
up and spin-down currents predicts the spin polarisation 
current to be proportional to the temperature gradient, 
corresponding to a spin Soret effect. 
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